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Abstract 

We are interested in the problem of characterizing the correlations that arise when 
performing local measurements on separate quantum systems. In a previous work 
[Phys. Rev. Lett. 98, 010401 (2007)], we introduced an infinite hierarchy of conditions 
necessarily satisfied by any set of quantum correlations. Each of these conditions could 
be tested using semidefinite programming. We present here new results concerning 
this hierarchy. We prove in particular that it is complete, in the sense that any set of 
correlations satisfying every condition in the hierarchy has a quantum representation in 
terms of commuting measurements. Although our tests are conceived to rule out non- 
quantum correlations, and can in principle certify that a set of correlations is quantum 
only in the asymptotic limit where all tests are satisfied, we show that in some cases it 
is possible to conclude that a given set of correlations is quantum after performing only 
a finite number of tests. We provide a criterion to detect when such a situation arises, 
and we explain how to reconstruct the quantum states and measurement operators 
reproducing the given correlations. Finally, we present several applications of our 
approach. We use it in particular to bound the quantum violation of various Bell 
inequalities. 



1 Introduction 



The main goal of Quantum Information Science (QIS) is to understand the possibilities 
and limitations of the quantum formalism for information processing and communication. 
Research in QIS is concerned on one hand with the design of new protocols exploiting the 
transmission and manipulation of information encoded in quantum states (see for instance 
PQ). On the other hand, it seeks to identify the constraints on information processing imposed 
by the quantum formalism. For instance, various information tasks, such as unconditionally 
secure bit commitment, have been shown to be impossible in a quantum world [2]. 
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Figure 1: Local measurements on a system shared by two observers viewed as a black-box 
process. Alice chooses a measurement input X and obtains a measurement output a G X. 
Similarly, Bob chooses an input Y and receives an output b &Y. The behavior of the system 
is characterized by the joint probabilities P(a, b). 



A standard scenario in QIS, and which serves as a primitive for more complex protocols, 
consists of two distant, non-communicating parties, conventionally called Alice and Bob, who 
share a quantum system in a joint state p. Each party makes a measurement on his share 
of the state and obtains a classical outcome. On a phenomenological level, we may describe 
the situation by saying that the two parties have access to a black box (see Figure 1). When 
Alice inputs a measurement X into the box, she gets as output a measurement outcome 
a G X; similarly, when Bob inputs a measurement Y, he receives an output b G Y. The 
behavior of the box is completely characterized by the joint detection probabilities P(a,b). 
From now on, we simply call a behavior the set P = {P(a, b)} of all such probabilities. 

Though in the above scenario the parties are separated and perform local measurements, 
their outcomes a and b may be non-trivially correlated, in particular if the initial quantum 
state p is entangled. These correlated data can be exploited for different tasks, such as 
communication complexity [3] or key distribution [3]. From the perspective of QIS, it is thus 
meaningful to characterize which outcome correlations can, or cannot, be produced by two 
non-communicating quantum observers. The main problem with which we are concerned in 
this paper is thus the following: given a behavior P, do there exist a quantum state p and 
local measurements X and Y reproducing the outcome probabilities described by PI Note 
that we do not impose here any constraints on the dimension of the system shared by Alice 
and Bob, as we are interested in the most general set of correlations that can be obtained 
with quantum resources. 

The special case of classical observers is relatively well understood. The correlations 
obtained in this scenario coincide with the ones that can be achieved with shared random- 
ness, or using another terminology, with those that are described by local-hidden variable 
models [5] . For a given number of possible measurement inputs and outputs, the set of local 
classical correlations forms a convex polytope whose vertices correspond to all the possible 
deterministic assignments of outputs to inputs. It thus follows that linear programming can 
be used to decide if a given behavior is reproducible by two local classical observers [HI El- The 
facets of the classical polytope, which form the boundary of the classical region, correspond 
to the well-known Bell inequalities [8]. 

Our understanding of the general case of quantum observers, with which we are concerned 
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Figure 2: Schematic representations of the space of joint distributions P(a, b) (for fixed and 
finite number of possible inputs and outputs). L denotes the set of correlations that admit a 
local model; it is a polytope and membership in L can be decided using linear programming 
[SI El- NL is the global set that contains all (in particular non-local) correlations; it is again 
a polytope. The region accessible to quantum mechanics is Q. The quantum set is not a 
polytope, i.e., it does not have a finite number of extreme points. As represented on the 
figure, Q contains L and is a proper subset of NL. See [9] or [10J for more details. 

here, is more rudimentary. The difficulty lies in the fact that we do not have a practical 
characterization of the set of quantum behaviors and that this set cannot be described by a 
finite number of extreme points (see Figure 2) [9j. 

Apart from the QIS motivation, the problem of characterizing the set of quantum behav- 
iors is also of relevance from a fundamental perspective. Indeed, while Quantum Mechanics 
has been so far confirmed by plenty of experiments, we cannot exclude that someday it will 
be disproved. If some experimental data were inconsistent with the quantum model that we 
have for the experiment, would that however necessarily imply the breakdown of the whole 
quantum formalism? How could one exclude that there is no other quantum model explain- 
ing the observed data? The problem then is to establish experimentally testable conditions 
that can rule out the whole quantum structure, in a similar fashion as Bell's inequalities 
do for locally causal model. In this context, one is again confronted with the problem of 
characterizing the constraints on correlations imposed by the quantum formalism. 

One of the first researchers to study the characterization of quantum correlations was 
Tsirelson in 1980 [TTJ. Tsirelson got several important results for the case of measure- 
ments with binary outcomes [9J; most notably he derived the maximal quantum violation 
of the Clauser-Horne-Shimony-Holt (CHSH) inequality [T2]. More recently, the problem has 
attracted the interest of several researchers working in QIS [13J. Among the latest contri- 
butions, we point out the work of Wehner [TJj, who showed that part of Tsirelson's findings 
could be implemented using a relatively new numerical tool called semidefinite programming 
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(SDP). A short introduction to this technique is given in Appendix A, more details can be 
found in [T5]. Apart from Wehner's paper, there are several other articles using SDP tech- 
niques to bound the set of quantum correlations [TBI HT] - Most of these results deal with the 
case of two-outcome measurements. 

In a recent work [18] . we introduced a hierarchy of SDP tests to check if a given behavior 
admits a quantum representation. Compared to previous constructions, our method is com- 
pletely general as it can be applied to any number of parties, measurements, and outcomes, 
and is independent of the dimension of the quantum systems. In this work, we explore 
further the approach introduced in [18]. 

The basic idea behind our method is presented in Section [3j Instead of directly searching 
for a quantum state and measurement operators reproducing a given behavior — a computa- 
tionally highly difficult task, if not impossible if the dimension of the system is unbounded - 
we consider instead a family of weaker conditions. Each of our conditions amounts to ver- 
ify the existence of a positive semidefinite matrix whose structure depends on the general 
algebraic properties satisfied by quantum states and measurement operators. If one of our 
conditions is not satisfied, we can immediately conclude that the given behavior is not quan- 
tum. In Section SI we show that our family of tests can be organized as an infinite hierarchy 
of increasingly stronger conditions. We prove that in the asymptotic limit, our hierarchy 
is complete in the sense that any behavior that satisfies all the conditions in the hierarchy 
necessarily has a quantum representation in terms of commuting measurements. We show 
further that in some cases it is possible to conclude that a behavior is quantum after a 
finite number of steps only. We provide a criterion to detect when such a situation arises, 
and we explain how to reconstruct in this case explicit quantum states and measurement 
operators reproducing the given behavior. Based on these latter results, we then show how 
a slight modification of our tests allow us to reduce the problem of deciding if a behavior 
has a quantum representation with quantum systems of finite dimension d to a rank min- 
imization problem. Unfortunately, and contrary to SDP, there are no efficient algorithms 
to solve rank-minimization problems. In Section we present several applications of our 
method. We use it in particular to put upper bounds on the quantum violation of various 
Bell inequalities. We conclude with a discussion and some open questions in Section [6j 

2 Definitions 

2.1 Measurement scenario 

We consider measurement scenarios as illustrated in Figure 1. We assume that outputs 
corresponding to different inputs are labeled in a distinct way. Each output, say a of Alice, 
is thus uniquely associated to a single input X(a). We denote by A the set of all outputs 
of Alice and by B the set of all outputs of Bob. The inputs of Alice may be viewed as 
disjoint subsets of A, and those of Bob as disjoints subsets of B. A measurement scenario 
is thus specified by a quadruple (A, B, X, y), where X and y are partitions of A and B, 
respectively. 
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The measurement scenarios that we consider in this paper always involve a finite number 
of inputs and outputs, i.e., A and B are finite sets. A behavior P thus consists of a finite 
set of \A\ x \B\ joint probabilities: P = {P(a,b) : a G A, b G B}. For instance, in the case 
where Alice and Bob have each a choice between s different inputs that each yield one out 
of d outputs, a behavior consists of s 2 x d 2 joint probabilities. Except otherwise mentioned, 
we assume in the remaining of the paper that a measurement scenario (A, B,X,y) and a 
behavior P have been specified. Our aim is to determine if P represents a possible quantum 
process. 

2.2 Quantum behaviors 

Definition 1. The behavior P is a quantum behavior if there exists a pure state in 
a Hilbert space 7i, a set of measurement operators {E a : a G A} for Alice, and a set of 
measurement operators {Eb : b G B} for Bob, such that for all a G A and b G B 



with the measurement operators E satisfying 

1. E\ = E a and E\ = Eb (hermiticity) 

2. E a E- a = 5 a aE a ifX(a) = X{a!) and E b E~ b = 5 b iE b ifY(b) = Y(b') (orthogonality) 

3. Yla^x^a = an d Ylb& Eb = (completeness) 
4- [E a ,E b ] =0 (commutativity) 

The set of all quantum behaviors will be denoted by Q. 

The first three properties are necessary to ensure that the operators E a and Eb are 
projectors and define proper measurements. The fourth property simply expresses the fact 
that Alice and Bob perform separated measurements on the global state \ip). 

Note that more generally we could have defined a quantum behavior in terms of a mixed 
state p rather than a pure one and in terms of general measurements, also known as Positive 
Operator Valued Measures (POVM) [T], rather than projective ones. But remark also that in 
our definition we put no restrictions at all on the dimension of the Hilbert space. Since any 
general measurement on a given Hilbert space can be viewed as a projective measurement 
on a larger Hilbert space, and any mixed state p can be viewed as a subsystem of a larger 
system in a pure state \ip) pQ, the above definition turns out to be completely general. 

Property 3 implies that the marginal probabilities P(a) = J2beY P( a > an d P(b) = 
Saex P( a i are well-defined and independent of what is measured on the other side (i.e., P 
satisfies the no-signalling constraints). This also implies that in the above definition there is 
some redundancy in the specification of the operators {E a : a G A} of Alice since any one of 
them can be written as the identity minus the other ones. To simplify further the definition 
above, select an output a x G X for each input X and introduce the reduced output sets 
X = {a : a G X,a ^ ax} and A = [J x X. Introduce analogous sets Y and B for Bob. The 
following definition is then equivalent to Definition [TJ 
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Definition 2. The behavior P is a quantum behavior if there exists a pure (normalized) 
state in a Hilbert space 7i, a set of measurement operators {E a : a G A} for Alice, and 
a set of measurement operators {E b : b G B} for Bob such that for all a G A and b G B 

P(a) = {^\E a \^) 
P(b) = (Vi£ 6 |V> 
P(a,b) = (4>\E a E b \i:) (2) 

with the measurement operators satisfying 

1. E\ = E a and E\ = E b (hermiticity) 

2. E a E- a = 5 a - a E a ifX(a) = X(a') and E b E b = 5 bb E b ifY(b) = Y(b') (orthogonality) 

3. [E a ,E b ] =0 (commutativity) 

It is clear that any behavior satisfying Definition[I]also satisfies Definition[2j The converse 
statement is also true. Indeed given sets of operators {E a : a G A} and {E b : b G B} 
satisfying Definition [2l define the missing operators E ax and E by through E ax = 1— ^2 a& x E a 
and E by = 1 — J2beY Eb- It is then easy to see that the now-complete sets {E a : a G A} 
and {E b : b G B} satisfy Definition [U 

Before concluding this subsection, note that when dealing with finite dimensional Hilbert 
spaces, one tends to associate a tensor product structure to separated measurements. This 
leads to another set Q' of quantum behaviors, possibly equivalent to Q, and defined as 
follows. 

Definition 3. The behavior P belongs to the set of quantum behaviors Q' if there exists 
a pure state \ip) in a composite Hilbert space ® TLq, a set of measurement operators 
{E a : a G A} for Alice, and a set of measurement operators {E b : b G B} for Bob, such 
that for all a G A and b G B 

P(a,b) = (ij\E a ®E b \tlj) J (3) 
with the measurement operators E satisfying 

1. E\ = E a and E\ = E b (hermiticity) 

2. E a E- a = 5 a - a E a ifX(a) = X{a') and E b E b = 5 bb E b ifY(b) = Y(b') (orthogonality) 

3. Ylia&x^a, = & a an d J2b&Y Eb = (completeness) 

Clearly, Q' CQ. However, it is an open question whether these two sets are equal. In the 
special case of finite dimensional Hilbert spaces, they turn out to be identical [19J. In this 
work, we adopt Definition [H or equivalently Definition [2j partly because it is much better 
tailored to the structure of our construction. We will come back to the commutation versus 
tensor product issue in section 16.11 
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2.3 Sets of operators and sequences 

In this subsection, we introduce a few other definitions that will be needed later on. 

Let £ denote the set of projectors appearing in Definition [TJ i.e., £ = {E a : a G A} U 
{E b : b G B}, and £ denote the set of projectors of Definition [2] plus the identity, i.e., 
£ = 1 U {E a : a G A} U {E b : 6g5}. 

Let (9 = {Oi, . . . , O n } be a set of n operators, where each Oi is a linear combination of 
products of projectors in £ . Thus O is a finite subset of the algebra generated by £. Note 
that we can equally well define the set O in terms of £, since £ and £ are equivalent up to 
linear combinations. Define T{0) as the set of all independent equalities of the form 

{F k )^\0\0^) =g k {P) k = 1, . . . ,m (4) 

which are satisfied by the operators Oi, where the coefficients gu{P) are linear functions of 
the probabilities P(a,b): 

9k(P) = (g k )o + J2(9k) ab P(a,b) (5) 

a,b 

and where is the state appearing in Definition [2j These equations are the ones that for- 
mally follow from the definition of the 0;'s, the relation (T5]), and properties 1-3 of Definition 
[2j Each set of operators O define such a collection of equations. As an example of equation 
of the form (TJJ, suppose that the set O contains the operators {Ok}f =1 = {E b E a S : a G X}, 
where S is some arbitrary operator in the algebra generated by £, and also contains the 
operator O d+1 = E b S. Then £Li 0\O k = J2 aeX (E b E a S^ E b E a S = £ aeX &E a E b E b E a S = 
J2 aex &E b E a E b S = &E b E b S = O j d+1 O d+1 , and thus ELi^l^O^) - ^\0 ] d+1 O d+1 \^) = 
0. 

Let a sequence S be a product of projectors in £. Examples of sequences are E a and 
E a E a iE b . Note that some sequences may correspond to the null operator, for instance, 
E a E' a = if a 7^ a', and X(a) = X(a'); in the following, when we speak of a sequence, we 
always mean a non-null sequence. The length \S\ of a sequence is the minimum number 
of projectors needed to generate it. For instance \E a E b E a \ = \E a E a E b \ = \E a E b \ = 2. By 
convention, the length of the identity operator is |1| = 0. We define S n to be the set of 
sequences of length smaller than or equal to n (excluding null sequences) . Thus 

5 = {1} 

51 = S U {E a : a G A} U {E b : be B} 

5 2 = S US 1 U {E a E a , : a, a G A} U {E h E v : b, b' G B} U {E a E b : a G A,b G B} 
Ss = ■ ■ ■ 

It is clear that Si C S 2 Q ■ ■ ., and that any operator Oi G O can be written as a linear 
combination of operators in S n for n sufficiently large. 
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3 Basic idea of our method 



The following proposition associates to each set of operators O satisfying Eqs. (J4j) a condition 
that restricts the possible correlations that can arise between two quantum observers. 

Proposition 4. Let O be a set of operators and T{0) the set of equations of the form (jlj) 
satisfied by operators in O. Then, a necessary condition for a behavior P to be quantum is 
that there exists a complex hermitian n x n positive semidefinite matrix T y whose entries 
Tij satisfy 

Y,{Fk) t] T l0 =g k {P) fc = l,...,m (6) 

ij 

Moreover, if the coefficients F k and g k in (jlj) are real, we can take V to be real as well. 

Proof. If P is quantum, there exist a state \ip) and projectors E a and Ej, as in Definition [21 
and therefore there also exist operators satisfying the relations (jlj). Then simply define 
the entries of the matrix F through 

T l3 = ^\0\0^) (7) 
Clearly, V satisfies ([6]). Moreover, it is positive semidefinite since for all v G C n 

v^Tv »;r, ; O = X>,*M0}0iM«i = WV^V\iP) > (8) 

ij ij 

where V = J2j v jOj- 

If the coefficients F k and g k in (jlj) are real, redefine V as (T + T*) /2. Then T still is 
positive semidefinite and satisfies ([6]). □ 

We will call a certificate associated to O to any nxn positive semidefinite matrix T satis- 
fying the linear constraints As an illustration, we now give two examples of application 
of Proposition |H 

Example 1 

Consider a measurement scenario where Alice has a choice between two measurements, X = 1 
or X = 2, to perform on her subsystem, and where both measurements yield binary outcomes 
with values ±X. Likewise, Bob has a choice between two measurements, Y = 3 or Y = 4, 
with outcomes ±Y. 

The single-party measurement averages Cx = P(+X) — P(—X) and Cy = P(+Y) — 
P(—Y) together with the two-party correlation functions Cxy = P(+X, +Y) +P(— A, —Y) 
— P(+X, —Y) — P(— X, +Y) fully determine the response of the joint system of Alice 
and Bob. The observed data are thus characterized by the eight numbers {C\, C 2 , C 3 , C 4 , 
012, C*i3, C23, C24} which are equivalent to the knowledge of the entire set of probabilities 
P(±X,±Y). 



S 



Criterion 5. If the data observed by Alice and Bob represent the response of a quantum 
system, there exists a real symmetric 5x5 positive semidefinite matrix T >z of the form 



r 




(9) 



1 v 



V 1 / 



where u andv are arbitrary entries. (We have only given the upper triangular part ofT since 
it is symmetric.) 

Proof. If the data observed by Alice and Bob represent the response of a quantum system, 
there exist a state two projectors E±x associated to each of the two measurements 
X = 1, 2 of Alice and two projectors E±y associated to each of the two measurements 
Y = 3,4 of Bob. Let O = {oo, <7i, . . . , 04} where 00 = 1 is the identity operator and 
(jj = E +i — E-i (i — 1, ... ,4). It is easily verified from Eqs. ([I]) and properties 1-4 that these 
operators satisfy the equalities 



which are the counterparts of Eqs. (jlj). It immediately follows that the associated 5x5 
matrix = {ip\a\a has the form Q). It can be taken real if we further redefine T as 



Example 2 

Consider the case where Alice and Bob have a choice between s different measurements 
that each yield one out of d possible outcomes. The s measurements of Alice are labeled 
X = 1, . . . , s and her m = s x d possible outcomes are labeled a — 1, ... ,m, where outcomes 
in the range 1 + (k — l)d, . . . ,kd belong to the measurement X — k. Analogously, the 
s measurements of Bob are labeled Y = s + 1, ... ,2s and his m = s x d outcomes are 
b = m + 1, . . . , 2m, where again outcomes in the range 1 + (k — l)d, . . . ,kd belong to the 
measurement Y = k. 

This measurement scenario is characterized by the m 2 joint probabilities P(a,b). 

Criterion 6. If the set of m 2 probabilities P(a,b) admits a quantum representation, there 
exists a 2m x 2m real symmetric positive semidefinite matrix T y of the form 




1 



0,...,4 



(10) 

(11) 
(12) 




(r + r*)/2. 



□ 




(13) 
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where the submatrix P is the m x m table of probabilities with entries P a b = P(a,b), and 
where the submatrices Q and R satisfy 



Qaa' = 5aaiP(0) if X{a) = X(d') (14) 

R w = 5 w P{b) if Y(b) = Y(b') (15) 

Proof. If the measurement scenario is a quantum measurement scenario, there exist a quan- 
tum state m projectors E a for Alice, and m projectors E b for Bob satisfying the prop- 
erties of Definition [TJ Consider the set O = £ = {E\, . . . , E m , E m+ %, . . . E 2m } consisting of 
the m operators of Alice and the m ones of Bob. They satisfy the equalities 

(V|#aW> =P(a,b) 

(iplEaKW) = Saa'P(a) if X(a) = X(a') 

{ip\E b E' b \iP) = 5 w P(b) if Y(b) = Y(b') (16) 

as implied by Eqs. ([T|) and property 2. It immediately follows that the certificate V associated 
to O has the form (USD- □ 



Note that the matrix ( 1131) can be thought of as a table of probabilities, where is the 
probability to obtain the two outcomes i,j G {1, . . . ,2m}. The only entries of this matrix 
which are not specified are the entries Y aa i = Q aa i associated to different measurements of 
Alice, X(a) ^ X(a'), and the entries = Rw associated to different measurements of Bob, 
Y(b) Y(b') . This is coherent with our interpretation of T since in a quantum scenario these 
entries correspond to non-commuting measurements performed on the same subsystem and 
are thus not jointly observable. Nonetheless, if the correlations P(a, b) have a quantum origin 
it is possible to assign a numerical value to these undetermined entries, namely (ip\E a E a > 
and {^\E b E v \^% such as the overall matrix (fl~3|) is positive semidefinite. 

3.1 Testing the existence of a certificate with SDP 

Checking the existence of a certificate T, such as the ones given in Examples 1 and 2, can 
be cast as a semidefinite program. Indeed it amounts to solve the following problem 

maximize A 

subject to tr (F%T) = g k (P) k = 1, ...,m (17) 

r - xi y o 



which after some elementary manipulations can be put in the form (1A-46j) . A positive 



solution A > to the above problem implies that there exists a positive semidefinite matrix 
r y XI y compatible with the linear constraints (jH]). A strictly negative solution A < 
means that any matrix F compatible with ([6]) is necessarily negative definite and thus that 
the given behavior P does not represent the outcome of a quantum experiment. 



1 Or the real part of these expressions, if we take T real. 
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As mentioned in the Appendix, there exist many available programs to solve problems 
of the type (TTTj) . Such programs solve these problems both in their primal and dual forms. 
The dual of flT7} is 

minimize ^VkgkjP) 

k 

subject to F(y) = J2ykFZhO (18) 

k 

J> fc tr(if) = l 

k 

If a program returns a negative solution for the primal for a given behavior P*, it also yields 
a dual feasible point y such that J2kVk9k(.P*) < 0. This dual feasible point provides a 
proof that the given behavior P* is not quantum; it can be interpreted as a quantum Bell 
inequality violated by P* in the sense that ^2 k ykgk(P) > is a linear inequality satisfied 
by all quantum probabilities. Indeed, the coefficients g k (P) defined in ([5]) depend linearly 
on the probabilities P(a, b), and thus the expression J2Vk9k(P) is a linear expression in the 
probabilities P(a, b). Moreover, from the second line of (fTSl . we deduce that for all behaviors 
P having a positive certificate r y 0, in particular, for all quantum behaviors, this linear 
expression is positive: Y^kVk9k(P) = Sfc 2/fctr (Fj T) = tr (F(y)T) > since r y 0. The 
behavior P*, however, violates this inequality, J2 k yk9k(P*) < 0, which demonstrates that it 
does not belong to Q. 

3.2 Equivalence between certificates 

Each set O of operators that we can write down yields a different condition satisfied by 
quantum theory. However, not all conditions built in this way are independent, as the 
following lemma shows. 

Lemma 7. Let O and O' be two sets of operators such that every operator in O' is a linear 
combination of operators in O. Then, the existence of a certificate V associated to O (for a 
given P) implies the existence of a certificate V associated to O' . 

Proof. By hypothesis, every operator 0\ E O' can be written as 0[ = J2 k C ik O^ where O k E 
O. Define then = ^2ki^kJ^kiCij. It is clear that T' satisfies the equalities (jSJ) associated 
to O', given that F satisfies the ones associated to O. We also have that T' = C^YC y 0, 
and thus V is a certificate associated to C. □ 

The criterion of Example 2 for s = 2 and d = 2, for instance, is equivalent to the one 
of Example 1, because the set of eight operators {E + %, ELi, E +2 , E_ 2 , E +3 , E_ 3 , E +4: , -EL4} is 
linearly equivalent to the set of five operators {o~ , 0"i, cr 2 , <7 3 , cr 4 }. 

In numerical implementations, we have of course always interest to use a criterion based 
on a set O of linearly independent operators so as to minimize the size of the matrices 
involved. Note also that to check systematically all the conditions that follow from our 
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approach, it is sufficient to check the ones associated with the sets S n defined in Subsection [23] 
since they generate by linear combinations all other possible operators. This point is made 
more precise in the next section. 

4 A Hierarchy of necessary conditions 

Motivated by the above lemma, define a certificate of order n, denoted T n , as a certificate 
associated to the set of operators S n . A certificate of order n is thus a \S n \ x \S n \ matrix 
and to index its row and columns we will use symbols that are in direct correspondence with 
the elements of <S n . Sequence operators S, E a , E a S, and 1 will be associated with row or 
column indices, s, a, as, and 1, respectively. We define the length \s\ of an index s to be the 
length \S\ of the corresponding sequence S. A certificate T n is thus a matrix with entries 
{T™ t : \s\, \t\ < n}, which according to the proof of Proposition H] may be interpreted as 
T™ t = (%l)\S^T\%l)) if P is a quantum behavior. 

From Proposition H] and the definition of the set S n , we deduce that T n is a real positive 
semidefinite matrix that satisfies the linear equalities 

r?,i = 1 , TL = P(a), Tl b = P(b), Tl h = P{a,b) (19) 

for all a G A and b G B, and 

r: }t = r: tV i£&T = rtv (r? it = o if5*T = o) (20) 

for all \s\, \t\, \u\, \v\ < n. Here the relations S*T = WV (or S^T = 0) are the ones that 
follow from properties 1-3 of Definition 2. For instance, T™ ba = r" afe , and r^ fea , = if 
X(a) - X(a'). 

As we mentioned earlier, 5i C «S 2 C . . . C «S n C . . and thus the family of certificates 
r 1 , T 2 , . . . , T n . . . , represents a hierarchy of conditions satisfied by quantum probabilities, 
where each condition in the hierarchy is stronger than the previous ones. Moreover, since 
in the limit n — > oo the linear span of S n coincides with the entire algebra of operators 
generated by S, this hierarchy embraces, according to Lemma [7J all the conditions that can 
be built from our approach. The strategy that we propose to verify the quantum origin of a 
given behavior P is thus the following. Check first if there exists a certificate T 1 of order 1 
associated to P. If there is no such certificate, we can conclude that the behavior P is not 
quantum, otherwise check the existence of a certificate T 2 of order 2. Repeat the procedure 
with certificates of increasing order as long as the behavior P satisfies the previous tests. 

A geometrical interpretation of our hierarchy is given in Figure [31 where Q n denotes the 
set of all behaviors P for which there exists a certificate of order n. 
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Figure 3: Geometrical interpretation of our hierarchy. Q is the set of quantum behaviors. 
Q n denotes the set of all behaviors for which there exists a certificate of order n. Testing 
the existence of a certificate of order n amounts to determine if a given behavior P belongs 
to Q n . Certificates of higher order provide a more accurate approximation of the quantum 
set Q, but are more demanding from a computational point of view. 

4.1 Sufficiency of the hierarchy 

We now show that our hierarchy is complete in the sense that lim n ^ 00 Q n = Q, or in other 
words, that any non-quantum behavior P necessarily fails one of our conditions at some step 
in the hierarchy. 

Theorem 8. Let P be a behavior such that there exists a certificate T n of order n for all 
n > 1. Then P belongs to Q. 

Proof. The proof proceeds in two steps. We first show that the sequence of certificates T n 
admits a proper limit lim^oo T n — > T 00 . We then construct from the matrix T°° a quantum 
state and quantum operators acting on a (possibly infinite-dimensional) Hilbert space TC 
that reproduce the behavior P. 

Note first, as shown in Appendix B, that all the entries {r™ t : \s\, \t\ < n} of the matrices 
T n are bounded by 1, i.e., \T™ t \ < 1. Now, complete each matrix T n with zeros to make it 
an infinite matrix r n with entries {T™ t : \s\, \t\ = 0, 1, . . .}; we can then view the matrices 
T n as infinite vectors in (the normed space of all bounded sequences u = (u±,U2, ■■■), with 
norm given by ||lt||cc = suPi As the sequence {f n : n = 1,2. . .} belongs to the unit 

ball of loo, it admits, by the Banach-Alaoglu theorem, a subsequence {n^} that converges 
in the weak-* topology to a limit T ni —>■ T 00 when % — > oo [20]. This implies in particular 
pointwise convergence, i.e., 

hm f$ - r~ (21) 

for all s, t. From the pointwise convergence, we deduce that T°° satisfies Eqs. f jT9l) and (1201) 
for all s, t, u, and v. Moreover, let denote the submatrix of T n corresponding to the 
entries {f J t : \s\, \t\ < N}. Since h for all n and N, (ETJ) implies that y for all 
iV=l,2,..'. 
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In the remaining of the proof, we construct from the matrix r°° a state \4>) and operators 
{E a : a G A} and {E b : b G B} satisfying the properties of Definition 2. 

The fact that >z for all N implies that there exists an infinite family of vectors 
{\v s ) : \s\ = 0, 1, 2, . . .} whose scalar products reproduce the entries of r°°, i.e., 

= (v s \v t ) (22) 

for any s, t of length \s\, \t\ = 0, 1, 2 . . . One way to establish this fact is through a sequential 
Cholesky decomposition of the matrices [21] . 

We now take as our Hilbert space 7i the vector space spanned by the vectors \v s ) and 
define, for all a G A, projectors E a as follows 

E a = proj (span{|u as ) : \as\ = 1, 2 . . .}) (23) 

where proj(V) is the projector on the subspace V. Since (E a S)^ E a >T = 5 a , a 'S^E a T when 
X(a) = X(a'), it follows from ( 120]) and ( |22i) that (u a s|^a't) = ^^'(^sl^at), which in turn 
implies that 

E*\Va>s) = 5aa>\v as ) if X(a) = X{a!) (24) 

An immediate consequence of this is that 

E a E a , = b aa ,E a ifX(a)=X(a') (25) 

i.e, that the operators {E a : a G X} form an orthogonal set of projectors. They thus satisfy 
properties 1 and 2 of Definition 2. 

Let us now examine the action of E a over an arbitrary vector \v s ). We find that 

E a \v s ) = E a \v as ) + E a (\v s ) - \v as )) 
= \v as ) + E a (\v s ) - \v as )) 

= \v as ) (26) 

The last identity follows from the fact that (v at \v s ) — (v at \v as ) = which can be deduced from 
(120]) , (1221). and the relation {E a T) ] S - {E a T) ] E a S = 0. Property fl2JD implies in particular 
that 

E a \vi) = \v a ) (27) 

By repeating the above construction, we can build operators {E^ : b G B} for Bob that 
satisfy properties analogous to ff25l) . ff26|) . and (T27j). From fT26l) . fT27j) . and the corresponding 
relations for Bob, we deduce by induction that 

S\ Vl ) = \v s ) (28) 

for any sequence S of length \S\ = 0, 1, 2 ... of the projectors {E a : a G A} and {Eb : 6 G B}. 
Combining Eqs. (J22|) and fT28|) . we find that 

r- = (0|stf |0) (29) 
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where we have defined \<f>) = \vi). Note that \<f>) is a normalized vector since {4>\4>) = = 1. 
Eqs. (1291) together with (1191) imply that the state |0) and the operators E a and E b satisfy 
Eqs. ©. 

It now remains to verify property 3, i.e., that [E a , Ej,] = 0. From the relation (E a Sy E b T '— 
(E b SyE a T = 0, the properties ([20]) satisfied by r°° and ([29]), we deduce that 

(^[E a ,E b ]f\ ( f ) ) = (30) 

for any sequences S, T of length \S\, \T\ — 0, 1, 2 . . . As the vectors S\<p) and T\<p) span the 
support of the operators E a and fl30|) implies that the commutator [S a , £y is equal to 
zero. □ 

Corollary 9. Q is a closed set. 

Proof. From Theorem [8] we know that Q = H^i Q % - As eacn °f the sets Q l is closed, its 
infinite intersection must be a closed set as well. □ 



4.2 Stopping criteria and extraction of quantum state and mea- 
surements 

Our hierarchy of conditions characterizes the quantum set Q in an asymptotic limit. Testing 
only a finite number of our conditions may at most allow us to conclude that a given behavior 
does not belong to Q (more precisely, testing the conditions up to the n th step in the hierarchy 
allows us to detect all behaviors that do not belong to Q n ). We now show that in certain 
cases, it is possible to conclude at a finite order n in the hierarchy that a given behavior P 
does belong to Q . In this case, we can also recover from the certificate T n the quantum state 
and the measurements E a and E b reproducing the behavior P. 

Let r n be certificate of order n associated to the behavior P. Fix a pair of inputs X, Y 
and consider the set of all sequences of the form E a E b S where a G X and 6 6 7, together 
with all sequences of length n — \. Let Jx,y be the set of indices associated to such sequences. 
Define as the submatrix of T n with entries {T™ t : s,t G Jx.y}- If 

rank(r^y) = rank(r n ), (31) 

for all X, Y, then we will say that the certificate T n has a rank loop. 

Theorem 10. A behavior P has a quantum representation of finite dimension if and only if 
P admits, for some finite N , a certificate T N of order N with a rank loop and rank(r Ar ) < d. 

Here, by a representation of dimension d, we mean that there exist a quantum state 
|\&) G TC and a set of operators {E a ,Eb G B(H)} satisfying the conditions of Definition [2] 
for some Hilbert space 7i of finite dimension dim(7Y) = d. We denote by Qd the set of all 
behaviors having a <i-dimensional quantum representation. 
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Proof. We first prove the if P has a finite dimensional representation, there exists a certificate 
of order n with a rank loop. As P G Qd-, there exist a state \<p) E Ti and projective 
measurements G B(TC), as in Definition 2, for some Hilbert space of dim(7-^) = d. The 
matrix T n with entries T™ t = (</>| ( S''T'|0) for all S, T G <S n is clearly a certificate of order 
n associated to P. Because S n C «S n+ i, T n is a submatrix of r n+1 for any n, and thus 
rank(r n ) < rank(r n+1 ). On the other hand, the space generated by the vectors S\<p), S 
being an arbitrary sequence has a dimension less or equal than dim(7i) = d and therefore 
rank(r n ) < d for all n. These two conditions imply that there exists an N such that 
rankfr^) = rank(r Af+1 ) < d. It follows that rank(r^+ 1 ) = rank(r 7V+1 ), for all X,Y, and 
thus that T N+1 has a rank loopJE 

Let us now prove the converse statement. Suppose thus that P admits a certificate T N 
with a rank loop and satisfying rank(r iV ) = d. Similarly to Section 14.11 we can perform a 
Cholesky decomposition of T N to write T^ t = (v s \vt) for some finite set of vectors {\v s ) : 
l s l < X}, whose span is a vector space of dimension at most d. Again as in Section fl~Tl we 
can then define a set of operators A = {E a : a G A} as 

E a = proj (span {\v as ) : \as\ < N}) . (32) 

It is easy to see that these projector operators satisfy (J2"5|) . and using the same arguments 
as in section 14.11 one can see that they also fulfill 

E a \v s ) = \v as ), (33) 

for | as | < N. In an analogous way, we build operators B = {E^ : b G B} for Bob with the 
same properties. It is then immediate that (vi\S^T\vi) = Tf t , for sequences j^l, \T\ < N. 
In particular, (vi\E a Eb\vi) = P(a,b). The operators in A and B thus satisfy Eq. (j2j) and 
conditions 1 and 2 of Definition El It remains to show that they also satisfy condition 3, 
i.e., commutativity. 

Take any quadruple a, b, X, Y such that a G X, b G Y, and consider the set of sequences 
S XY = {S : \S\ < N - 1 or S = E a E b S', with a G X, b G Y, \S\ < N}. Then, for any pair 
of sequences S, T G Sxy, 

(v s \E a E b - E b E a \v t ) = Y N asM - Tl M = 0, (34) 

where the last equality comes from the constraints imposed on the certificate T N by the 
operator identity E a E b T — E b E a T = 0. On the other hand, from condition (I3T1) . we 
have that 

span({|^> : \s\ < N}) = span({|^) : S G S XY }). (35) 

Since the first set of vectors spans the support of the operators E a , E b , relation fl34l implies 
that E a , E b commute. As this holds for any quadruple a, 6, X, Y, it follows that P G Qd- □ 

2 What we have proven is that P has a, in general, complex rank looped certificate. To see that P also 
has a real rank looped certificate, note that, for any set n, Re(r n ) is also a valid certificate for P. On the 
other hand, rank(Re(r™))< 2 • dim(7i), so the previous arguments can be applied to Re(r n ). 
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Corollary 11. Let P be a behavior corresponding to a bipartite system where Alice's (Bob's) 
measurements have oIa (ds) possible outcomes and such that each of the probabilities satisfies 
P(a, b) > 0. Let T 2 be a certificate of order 2 compatible with this behavior. Then, rank(T 2 ) = 
dAds implies that P G Qd, with d = dAds- 

Proof. If P(a,b) > 0,Va G A, V6 G B, then, for any pair of measurements X,Y, the dAds 
vectors {\v s ) : s — 1, E a , Ef,, E a Eb : a G X,b G Y} can be shown to be linearly independent. 
This, together with the fact that the rank of the whole matrix is equal to d^ds, implies that 
r 2 has a rank loop. □ 

The above theorem says that if our SDP outputs a certificate T with a rank loop, we 
know that P belongs to Qd, with d = rank(r). Moreover, from the proof of Theorem [TU1 it 
is not difficult to see that we can even reconstruct the state and measurements E a and 
Eb that yield this finite-dimensional representation. 

Given a behavior P admitting a quantum representation of dimension d, there may be, 
however, different certificates of order n compatible with P, including some without rank 
loops. We have no guarantee that our SDP will output a certificate that has a rank loop, 
and thus in general we cannot guarantee that our hierachy of SDP tests will stop after a 
finite number of iterations. 

In view of this, it would be useful to incorporate some rank minimization techniques in 
the implementation of our hierarchy. That is, when checking the existence of certificates of 
order n, we would like as well to minimize the rank of the corresponding matrices. Indeed, let 
T™ be the certificate of order n for P with minimum rank, and consider the series T 1 , T 2 , T 3 , ... 
If rank(f™ +1 ) ^ rank(f n ), then rank(f n+1 ) > rank(f n ) + 1. This, together with the fact 
that rank(r n ) < d for all n, implies that there exists some N < d such that rank(r iV+1 ) = 
rank(r JV ). On the other hand, for all X, Y 

rank(f w ) < rank(f^+ 1 ) < rank(f JV+1 ), (36) 

and so rank(f 7V+1 ) = rank(f^y 1 ), i.e., f N+1 has a rank loop. 

Unfortunately, there are no known efficient methods to solve rank minimization of pos- 
itive semidefinite matrices with linear constraints. There are, however, heuristics [22] that 
typically arrive at the optimal solution in just a few iterations. 

5 Applications 

In this section, we present several applications of our method. We first derive simple analytic 
conditions that are satisfied by all quantum probabilities involving two measurements with 
two possible outcomes. We then show how to apply our method to establish upper bounds 
on the quantum violation of Bell inequalities. 

From a general perspective, the hierarchy of necessary conditions that we have introduced 
represents a systematic way of getting better and better approximations to the set of quantum 
correlations. Moreover, these approximations are nicely characterized in terms of semidefinite 
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constraints. Our method can thus be useful in any kind of optimization problem over this 
set. This is particularly true when we want to optimize the violation of Bell inequalities 
since they are linear functions of the behaviors and thus the entire optimization problem can 
be cast as a SDP. 

Although the applications that we present here are restricted to a bipartite scenario, our 
method can also be applied to a multipartite scenario, e.g. see [23] . 



5.1 Analytic conditions for quantum behaviors with two inputs 
and two outputs 

Consider the measurement scenario described in Example 1 of Section [3l involving two mea- 
surements with two possible outcomes for each observer. As we showed, a necessary condition 
for a behavior to be quantum in this scenario is the existence of a positive semidefmite matrix 
of the form ([9]). This condition corresponds to the first one in our hierarchy and thus char- 
acterize the set of behaviors Q 1 . In the following, we provide an analytic characterization of 
this set. The conditions that we obtain can be interpreted as the quantum analogues of Bell 
inequalities. 

We make use of the following two lemmas: 

Lemma 12. (Schur's lemma) 121^ Let M be a matrix such that 

M=( ?Uo, (37) 



Q R . 

with PyO. Then, M y if and only if R - Q T P- 1 Q h 0. 
Lemma 13. Let M. z % be a real symmetric matrix of the form 



( 1 Z X\ x 2 \ 

1 x 3 x 4 

1 t 

V 1 I 



(38) 



with \xi\ < 1, i — 1,2,3,4. Let f(x\, x 2 , £3, £4) = arcsin(xi) + arcsin^) + arcsin(x 3 ) — 
arcsin(a;4). Then, there exists a pair of values (z,t) such that M z>t >z if and only if 

\f(x 1 ,x 2 ,x 3 ,x i )\ < 7T (39) 

for all possible permutations of x%, x 2 , £3, x±. 

Proof. See Appendix E. □ 

Now, apply Schur's lemma to matrix ([9]), taking the upper block to be P = 1 >- 0. It 
then follows that the positivity of ([9]) is equivalent to the positivity of the matrix V given 
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by 



V 



(I - C\ 



\ 



u 



C1C2 



C13 — C1C3 

i-cl 



(7x4 

C24 
t> - 



- c 1 c 4 \ 

- C2C4 
C3C4 



(40) 



/ 



Note that we can restrict our analysis to the case where all the elements in the diagonal 
are strictly positive. Indeed, if a diagonal element is equal to zero, the corresponding mea- 
surement, say by Alice, is deterministic, i.e., it always returns the same outcome. Then, 
Alice is left with one effective measurement (at most) and there always exists a classical, 
hence a quantum, model for this type of scenario. Suppose thus that all the diagonal el- 
ements of r' are different from zero. Multiplying V on both sides by the diagonal matrix 



(i-^r 1/2 ( 



,4), we obtain a matrix of the same form as the one of Lemma 



M u = 

(fl3|) . Applying this lemma, we conclude, together with the previous observation, that a 
necessary and sufficient condition for a behavior to belong to Q 1 is either that there exist an 
% such that C? = 1 or that 



E 



arcsm 




CiCj 



i-c?)(i-c?; 



2 arcsin 



Cm — CkCi 



V(i-C2)(i-C7) 



< 7T 



(41) 



for all k = 1,2 and I = 3,4. This condition is of course only a necessary condition for 
quantum behaviors. 

Note that a weaker necessary condition for a behavior to be quantum follows from the 
positivity of 

/ 1 u C 13 C 1A \ 

1 C23 C24 
1 V 

\ 1 / 

which is simply a sub matrix of (Q. A direct application of Lemma [TBI implies that a behavior 
is quantum if | J^. arcsin(Cjj) — 2 arcsin (Cm)\ < for all k — 1, 2, i = 3, 4. This condition, 
which, as we said, is weaker than (jUJ), had previously been obtained in 0, [2H 123] . 



(42) 



5.2 Quantum violation of Bell inequalities 

Bell inequalities are constraints satisfied by all behaviors that originate from classical non- 
communicating observers. As mentioned in the Introduction, for a finite number of measure- 
ments and outcomes, the set of behaviors achievable using classically correlated instructions 
(shared randomness) defines a polytope, that is, a convex set with a finite number of extreme 
points (see also Figure [2]). It can then alternatively be completely characterized by a finite 
number of facets, which correspond to the well-known Bell inequalities [5]. A given behavior 
P thus admits a local classical model if and only if it satisfies all the Bell inequalities. In 
the space of behaviors, a Bell inequality can be viewed as a hyperplane that separates the 
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space in two regions. A generic Bell inequality can thus be written as 



I(P) = J2c ab P(a,b)<I ( 



c, 



(43) 



a.b 



where c a b are the real coefficients defining the inequality and Ic is the maximal value achiev- 
able by local classical points (and in particular which is attained by the extreme points lying 
on the facet defined by the Bell inequality). 

Since the work of Bell [8], we known that some quantum behaviors are incompatible 
with a local classical description, that is, that they violate a Bell inequality. This fact is 
often referred to as quantum non-locality. In spite of many years of work on quantum non- 
locality, there are no methods able to provide the maximal quantum violation of a general 
Bell inequality, or just non-trivial upper bounds to it An important exception already 
mentioned in the introduction is the (tight) bound derived by Tsirelson on the maximal 
violation of the CHSH inequality. 

Since our hierarchy of necessary conditions provides better and better approximations to 
the set of quantum correlations, it can be used to derive better and better upper bounds to 
the quantum violation of a Bell inequality. Actually, our proof of completeness guarantees the 
convergence to the maximal quantum value, that we denote by Iq. That is, by maximizing 
the value I(P) of a Bell inequality over the behaviors P(a, b) e Q n admitting a certificate of 
order n, one gets an upper bound I n to Iq. Clearly, we have that I± > I2 ■ • ■ > I n > • • • > Iq 
and lim^oo J n — > Iq. Even while we are only able to prove convergence to the quantum 
value in the asymptotic limit, the quantum value or a very good upper-bound to it can often 
already be obtained for a small relaxation order n, as we show in the following. 

The fact that Bell inequalities are linear functions of the joint probabilities P(a, b) signifies 
that we can cast the computation of these upper bounds as SDP. Indeed, note that for any 
certificate T n , we can write the value I(P) of a Bell inequality as I(P) = tr(/3 n r n ), where (3 n 
is a matrix whose elements are all zero but the entries corresponding to r™ 6 . For instance, 
in the case n — 1 one has (see ffT3"|) ). 



where C is the matrix whose elements are the coefficients c a b in ( )43l) defining the Bell in- 
equality. Therefore the calculation of I n amounts to solve the following SDP 



In the remaining of this subsection we illustrate this approach by applying it to several Bell 
inequalities. 

But before presenting these results, let us make two technical remarks. First, note that in 
the above optimization problems, we can in general consider certificates that are intermediate 




(44) 



maximize tr(/3„r n ) 

subject to tr (f% T n ) = g k (P) k = 1, m 
T n y 



(45) 
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between, say, a certificate of order 1 and a certificate of order 2. Such a certificate would 
be associated to a set of sequences of operators S satisfying Si C S C «S 2 . For instance, 
we could consider the set Si + ab = $i U {E a Eb, ■ a 6 i,i 6 5} consisting of Si together 
with all products of one operator of Alice and one for Bob (while S2 also contains product 
of two operators of Alice and product of two operators of Bob). The corresponding bound 
h+AB would then satisfy Ii < Ii+ab < ^2- In some cases this bound might already be useful 
while requiring less computational resources than J 2 . In the following, we will therefore 
also consider such bounds based on intermediate certificates. The notation that we use is 
obvious, for instance I 1+ ab+aa' b is the bound associated with the set of sequence operators 
Si+ab+aa'b = «Si U {E a E b } U {E a E' a E b }. Note that the rank loop conditions derived in 
Subsection 4.2 generalize to the case of intermediate certificates, see Appendix C. 

The second technical remark is that, as shown in Appendix D, the probabilities P(a, b) = 
corresponding to a certificate T n are guaranteed to be positive only for certificates of 
order n > 2 (or more generally for certificates associated with set of operators S ~D Si + ab)- 
Thus, when we maximize, as in (145)) . a Bell inequality over all behaviors for which there exists 
a certificate T 1 of order 1, it may happens that the bound I\ that we obtain correspond to 
a solution with negative probabilities. By explicitly adding to the SDP (14"5l) . the constraints 
Tafe — that probabilities must be positive^), we thus strengthen the upper bound Ii. In the 
remaining of this section, when we mention an upper bound obtained from a certificate of 
order 1, we always refer to this strengthened version. 

We start by analyzing the Collins-Gisin-Linden-Massar-Popescu (CGLMP) family of Bell 
inequalities introduced in [27] . These inequalities are defined in a bipartite scenario where 
the two observers can each make two measurements of d outcomes. We refer the reader 
to the original reference for the detailed description of these inequalities. The inequality 
corresponding to the case d — 2 is the CHSH inequality. The best known lower bounds 
on the quantum violation of these inequalities for d < 8 are those given in Ref. [28]. The 
upper-bounds that we obtained using our method are given in Table 1. Note first, that 
in the case d = 2 (CHSH) the first certificate already provides the actual quantum value, 
which is equal to the Tsirelson bound. For d larger than 2, the quantum value is recovered 
at the successive step corresponding to the certificate Y 1+AB . This can be seen by noting 
that the upper-bounds Ii+ab are equal to the lower bounds given in [28J. Alternatively, 
one reaches the same conclusion by noting that the stopping critera based on rank loops 
presented in Section 14.21 are satisfied. Thus, Y 1+AB , and therefore T 2 , is already enough 
to get the maximal quantum violation of CGLMP inequalities (at least until d = 8) and 
certificates T n with n > 2 are redundant. 

We have also considered other, perhaps less standard, Bell inequalities, like the one 
presented in [29] (see also [30]) for the case in which Alice performs two measurements, 
one of two outcomes and one of three outcomes, while Bob performs three two-outcome 
measurements. The results are summarized in Table 2. One can also get numerical lower 
bounds for the maximal quantum violation for fixed dimension. In the case of qutrits, the 
derived quantum violation is equal to 0.2532 [291 [30]. This is precisely the same value 

3 Adding such constraints leaves the optimization problem in a SDP form 
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d 


h 


h+AB 


Value 


Matrix Size 


Rank Loop 


Value 


Matrix Size 


Rank Loop 


2 


2.8284 


5 


N/A 


2.8284 


9 


Yes 


3 


3.1547 


9 


N/A 


2.9149 


25 


Yes 


4 


3.2126 


13 


N/A 


2.9727 


49 


Yes 


5 


3.2997 


17 


N/A 


3.0157 


81 


Yes 


6 


3.3378 


21 


N/A 


3.0497 


121 


Yes 


7 


3.3843 


25 


N/A 


3.0776 


169 


Yes 


8 


3.4115 


29 


N/A 


3.1013 


225 


Yes 



Table 1: Upper bounds on the violation of the CGLMP inequality derived from our con- 
struction. The local bound is equal to 2. The upper bound h+AB is already equal, up 
to numerical precision, to the lower bounds given in [28]. We also provide the size of the 
certificates in each case. Note that the rank loop conditions defined in Subsection (14. 2p are 
not applicable to certificates of order 1. 



Upper bound 


Value 


Matrix Size 


Rank Loop 


h 


0.3333 


7 


N/A 


I\+AB 


0.2653 


16 


No 


h+AB+AA'B 


0.2532 


22 


Yes 



Table 2: Upper bounds on the violation of the Is inequality derived from our construction. 
The local bound is equal to 0. The upper bound I\+ab+aa>b is already equal, up to numerical 
precision, to the lower bound obtained numerically for qutrits. We also provide the size of 
the certificates in each case. 



obtained when checking the last certificate of Table 2. This certificate then, or equivalently 
T 3 , already provides a tight bound on the maximal quantum violation. The same conclusion 
follows again by studying the rank of the matrices appearing in these certificates. 

Finally, we also applied our techniques to the Froissard inequality, also referred to as hz22 
inequality, given in [31, 32J. Again, we refer the interested reader to these references for the 
explicit form of the inequality. The best known quantum violation of this inequality is equal 
to 0.25 in the case of qubit systems, while the classical value is equal to zero. By applying 
our hierarchy of conditions to this inequality, one gets the upper bounds given in Table 3. 
Note that the values derived for T 2 and T 3 are quite close and that no rank loop is observed]^ 
It is remarkable that none of our upper bounds coincides with the best known lower bounds 
on the quantum violation, although they are very close to it. This may be because in the 
case of this inequality our hierarchy approaches more slowly the quantum solution, assuming 
it to be equal to 0.25. However, one cannot exclude that the maximal quantum violation of 

4 Rank loops should be considered in a cautious way. Indeed it is sometimes difficult to numerically 
distinguish a zero from a small eigenvalue. 
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Upper bounds 


Value 


Matrix Size 


Rank Loop 


h 


0.3333 


7 


N/A 


h+AB 


0.2515 


16 


No 


h 


0.25091 


28 


No 


h 


0.25089 


88 


No 



Table 3: Upper bounds on the violation of the 13322 inequality derived from our construction. 
The local bound is equal to 0. Interestingly, none of our tests coincides with the best known 
lower bound on the quantum violation, obtained for qubit systems. We also provide the size 
of the certificates in each case. 



this inequality is obtained for systems of dimension larger than two. Indeed, the existence of 
this type of inequalities has recently been proven in [33], EU [30] . Thus, a quantum violation 
close to 0.2509 is perhaps attainable beyond qubits. 

6 Discussion and open questions 

Characterizing the correlations attainable by quantum means is a fundamental problem in 
Quantum Information Science and, more generally, in Quantum Mechanics. To our knowl- 
edge, our construction represents the only available tool to tackle this problem with full 
generality: it applies to any number of parties, measurements and outcomes. Moreover, 
the first steps in our hierarchy are easily computable since they correspond to semidefinite 
programs of reasonable size. Our construction provides a systematic way of getting better 
and better approximations to the set of quantum correlations and can be applied, for in- 
stance, to identify correlations that do not admit a quantum representation, or to estimate 
the maximum quantum violation of Bell inequalities. 

In this work, after having presented in detail the hierarchy of necessary conditions already 
introduced in [18], we have (i) proven the completeness of the hierarchy, (ii) introduced a 
criterion based on ranks loops that can guarantee at a finite order in the hierarchy that a 
set of joint probabilities is quantum, and we have shown in this case how to reconstruct 
the quantum state and measurements reproducing these probabilities, (iii) presented several 
examples illustrating the usefulness of the method. Although our results are described in 
the bipartite case, they can easily be extended to the multipartite scenario. To conclude 
this work, we would like to go back to the commutation vs tensor product issue briefly 
mentioned in Subsection I2.2[ discuss the computational complexity of our approach, and 
then present several open questions related to the set of quantum correlations achievable 
with finite dimensional Hilbert spaces. 

We mention that it is possible to generalize the hierarchy presented in this work and 
our proof of convergence to other polynomial optimization problems with non-commutative 
variables [33]. A similar generalization was also recently introduced in J3B] to put upper- 
bounds on the entangled value of quantum multi-player games. We also mention that an 
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alternative proof of convergence of our hierarchy is possible using a result of Helton and 
McCullough [H7], as noted in [SHIES]. 

6.1 Commutation vs tensor product 

There are two possible ways to impose that Alice and Bob perform measurements on sepa- 
rated systems: through the condition that their measurement operators commute, or through 
a tensor product splitting of the whole Hilbert space. The two sets of quantum correlations Q 
and Q' associated with each possibility are defined in Subsection 12.21 Clearly, measurements 
that have a product form commute with each other, and thus Q' C Q. In the special case 
of finite-dimensional systems, one can in fact show that both definitions are equivalent, i.e., 
Q — Q' [T9l 139] (see also [10] ) . For infinite-dimensional systems whether they are equivalent 
or not is, at the moment of writing, an open question [19J. It is not even known if Q' is dense 
in Q. (Note that the statement in [9] that the two sets are equivalent is actually unproven 

ES]). 

One can debate which definition should be regarded as the proper one. Arguments in favor 
of the tensor product structure are presented in [19]. Here, we have chosen commutativity 
as this choice is consistent with the ethos adopted in this work. Indeed, our main objective 
is to characterize the set of correlations compatible with the general structure of Quantum 
Theory but imposing as few additional constraints as possible: we impose no restrictions on 
states, measurements, or even on the Hilbert space dimension. In this spirit it should then 
be pointed out that there exist in quantum field theories algebras of local (in the sense of 
commuting) observables that cannot be split in a tensor product structure^), yet in which it is 
possible to investigate the correlations that can arise between two separated observers, and in 
particular to study the amount by which Bell inequalities are violated [H] . By investigating 
the structure of the set Q defined through commutativity, we are sure to include also these 
examples and thus to deal with the most general correlations compatible with Quantum 
Theory. 

Of course, making the above distinction is only meaningful if Q and Q' happen to be 
distinct. But note that actually most of the results of this work are independent of the defi- 
nition chosen. As Q' C Q, all the necessary conditions satisfied by points in Q, in particular 
all the ones constituting the hierarchy, are also valid for Q'. The stopping criteria presented 
in Subsection 14.21 are associated with correlations achievable with finite-dimensional spaces, 
for which we known that Q = Q', and thus also apply to both cases. The unique distinction 
arises when we consider the asymptotic behavior of our hierarchy: as our proof of conver- 
gence to Q explicitly use infinite-dimensional systems, the hierarchy will also converge to 
Q' only if Q — Q' in the most general setting. But for all practical applications of our 
method where only a finite number of steps of the hierarchy are involved, in particular for 
all numerical applications, one choice of definition or the other does not make any difference. 
For instance, all the results presented in Section \5\ apply equally well to both cases. 

5 This result does not directly imply that Q ^ Q' since it could happen that all the correlations obtained 
by performing commuting measurements on states belonging to these spaces can also be realized in spaces 
with a tensor product structure. 
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Note that it is not surprising that the limit of the hierarchy tends to Q rather than Q', 
as the space separation between Alice and Bob's measurements appears in the hierarchy 
only in the form of constraints associated to the commutativity of these local observables. 
For example, since E a Et,E a i = E a E a >Eb, we impose that r™ ba , = r™ a , 6 for all n. If we insist 
that the hierarchy should tend to Q' rather than Q, it will probably be necessary to add 
new constraints associated to the tensor product structure. These constraints will have to 
reflect the (at the moment unproven) differences, at the level of operator algebra, between 
the commutation and tensor product case. 

6.2 Complexity of the hierarchy 

The computational complexity of our tests scales badly with the order n of the relaxation. 
For instance, in a measurement scenario with s inputs and d outputs, it is not difficult 
to see that the size of a certificate of order n is roughly (ds) n . The algorithms used to 
solve the semidefinite programs associated with such certificates have a running time that is 
polynomial in the size of the matrix defining the SDP. Thus, using semidefinite programming 
to decide if a certificate of order n exists requires a time exponential in n. 

Note, however, that the numerical results presented in Section 15.21 suggest that it might 
be sufficient, at least for some families of measurement scenarios, to consider relaxations only 
up to a bounded value n to characterize, or obtain an already good approximation, of the 
quantum set. Indeed, in the examples that we considered, when maximizing the violation 
of Bell inequalities we hit the quantum value, or obtained a very good upper-bound on it, 
already at the second or third step in the hierarchy. The suggestion that a finite number of 
steps of the hierarchy might already characterize, or approximate well, the quantum region 
turns out to be true in some particular case. For instance, for measurement scenarios with two 
outputs, a result of Tsirelson [11] implies that deciding if a set of correlators (i.e., a quantity 
such as the defined in Example 1 of Section 3) is quantum can exactly be decided through 
semidefinite programming, as noted by Wehner [2]. The semidefinite program considered 
by Wehner is a weaker version of the first step of our hierarchy. In [12], the authors show how 
for a certain family of measurement scenario, corresponding to unique games, the quantum 
set can well be approximated through semidefinite programming. The semidefinite programs 
considered in [12] correspond again to the first step of our hierarchjj^. 

If all these results suggest that our construction might indeed provide an efficient char- 
acterization of the quantum set for some particular quantum scenarios, we do not expect 
this to be true in full generality, as it has recently been shown, at least in the tripartite case, 
that calculating the maximal quantum violation of a Bell inequality is an NP-hard problem 

6 With the additional constraint, when maximizing the violation of Bell inequality, that the probabilities 
must be positive, as mentioned in Section 5.2. 
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6.3 Finite dimensional quantum systems 



In this work, we were mainly interested in characterizing the set of quantum behaviors 
without any bound on the dimension of the Hilbert space. We now present several open 
questions linked to the finite-dimensional case. 

• Consider all possible quantum behaviors of d outcomes where the number of measure- 
ments is arbitrary. Gill recently asked whether these correlations are attainable by 
measuring d-dimensional quantum systems [S]. The answer to this question is no, 
as shown in |33j for the case of three observers and in [301 EI! for bipartite systems. 
Actually, no finite dimension is sufficient to generate the whole set of quantum correla- 
tions of d outcomes for three parties, while the same result seems very plausible in the 
bipartite case [301 EI] • Consider however a scenario where the number of measurements 
is also finite. Are now all quantum correlations (exactly) attainable by measuring a 
finite dimensional quantum system? 

• Consider a measurement scenario with a finite number of inputs and outputs. It is 
easy to see that in the tensor product scenario discussed in Section 6.1. a quantum 
behavior can be approximated arbitrarily well using finite-dimensional Hilbert spaces 
(see for instance |36j). Does the same result hold in the commutative case? If yes, then 
combining this result with the fact that Q' = Q for finite-dimensional systems, would 
imply that Q' is dense in Q, and thus that our hierarchy converges to the quantum set 
Q' defined through the tensor product structure. 

• What is the structure of the set of quantum behaviors corresponding to a Hilbert space 
of fixed dimension dl Very little is known in this case, we even do not known if the 
corresponding quantum set is convex. 

• In relation with the above question, it would be interesting to understand how to in- 
corporate in our hierarchy a bound on the Hilbert space dimension. It is in principle 
always possible to decide if a behavior can be represented with a Hilbert space of given 
dimension through semidefinite programming [30J using known techniques of polyno- 
mial optimization H6j. The corresponding SDPs, however, are very demanding 
from a computational point of view, much more than the one obtained from our con- 
struction where we do not bound the dimension. Can one modify our construction to 
design more efficient methods to approximate the set of correlations corresponding to 
c?-dimensional quantum systems? 

As suggested by the results of Section 4.2, a possibility would be to incorporate a 
bound on the rank of our certificates. There are, however, to our knowledge no efficient 
methods to solve SDPs with rank constraints. Is there any efficient way to relax these 
rank constraints to obtain good approximations to the set of quantum correlations 
with finite dimension? 
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Appendices 

A Basics of semidefinite programming 

Semidefinite programming [15j is a subfield of convex optimization concerned with the fol- 
lowing optimization problem, known as the primal problem 

maximize tr(GZ) 

subject to trFiZ — C{ i = 1, ...,p (A-46) 

z y o 

The problem variable is the n x n matrix Z and the problem parameters are the n x n 
matrices G, Fi and the scalars q. A matrix Z is said to be primal feasible if it satisfies the 
conditions expressed in flA-46j) . 

For each primal problem there is an associated dual problem, which is a minimization 
problem of the form 

minimize c T x 

p 

subject to F(x) = ^ x^ - G h (A-47) 

i=l 

where the variable is the vector x with p components X{. The dual problem is also a semidef- 
inite program, i.e., it can be put in the same form as (1A-46I) . A vector x is said to be dual 
feasible when F(x) > 0. 

The key property of the dual program is that it yields bounds on the optimal value of 
the primal program. To see this, take a dual feasible point x and a primal feasible point 
Z. Then c T x - tr(GZ) = £Xi to(ZFi) x i ~ tr(GZ) = tr(ZF(x)) > 0. This proves that the 
optimal primal value p* and the optimal dual value d* satisfy d* < p*. In fact, it usually 
happens that d* = p*. A sufficient condition for this to hold is that there exists a strict 
feasible point of the primal problem [T5], that is, that there exists a matrix Z y that is 
primal feasible. Such a situation appears in the SDP problem ffT7|) . as for any matrix V 
satisfying the corresponding linear constraints, we can always take A small enough so that 
r-Al yO. 
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There exist many available numerical packages to solve SDPs, for instance for Matlab, 
the toolboxes SeDuMi [47J and YALMIP pEg] . Such algorithms usually solve both the primal 
and the dual at the same time and thus yields bounds on the accuracy of the obtained 
solution. 

B Certificates have bounded entries 

Proposition 14. Let T n be a certificate of order n for a behavior P. Then, \T™ t \ < 1, for 
all sequences S, T. That is, the set of all certificates of order n for P is bounded. 

Proof. Because T n >z 0, it just suffices to prove that all diagonal elements are smaller or 
equal than 1. Consider thus any 2x2 submatrix of r": 

(A-48) 

This submatrix must be positive semidefinite or, equivalently, its coefficients have to satisfy 
r™,r« > and r™ • I™ > r™ • I™. Now, take T = E a S. From the operator relation 
gt T = rptg = T t T; it foH ows that = = r&. This, together with the positivity 
conditions, implies that 

VI < T n ss , for T = E a S,V\S\ <n,aeA. (A-49) 

In particular, 

r^<r?i = i fOTalla - ( A - 5 °) 

And, obviously, the same relations hold replacing a's by 6's. By induction, it is straightfor- 
ward that < 1,VS, and, therefore, |r™ | < 1,VS,T. □ 

C Rank loop conditions for intermediate certificates 

We state here some results about rank loop conditions similar to the ones introduced in 
Subsection 14.21 and which hold for "intermediate certificates" such as those that we used in 
Subsection 15.21 to maximize the violation of Bell inequalities. 

Let us first define more precisely the certificates that we are considering here. Given 
a pair of measurements X, Y, denote by S n+ xy the set of sequences S n U {S G <S n+ i : 
S = E a E b S',a G X, b G Y}, i.e., S n+ xy is the set that contains all sequences of length 
n together with all the sequences of length n + 1 that are of the form E a EbS' for some 
a G X ,b G Y . It is thus intermediate between the set of sequences of length n and n + 1, 
as S n C S n+ xY Q <Sn+i- Given a vector n of positive integers n XY , define S u+ ab as the 
union of all sets S nxY+ xy- By abuse of notation, when n is an integer we interpret it as the 
vector (n,n, . . . ,n). With the notation that we have just defined, we have for instance that 
<Si+ab — {1, E a , E b , E a E b : a,b G A, B}, which is one of the set of sequences that we used in 
the numerical applications presented in Subsection 15.21 
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Given an arbitrary certificate T associated to a set of operators S and a vector n of 
positive integers such that S n+ xy Q S, denote by T n+X y the submatrix of T corresponding 
to the set of sequences S n+ xy- Define similarly T h+ ab- If there exists a vector N such that 

rank(r JVx¥+ xy) = rank(rjv+AB), (A-51) 

for all X, Y, then we will say that the certificate T has a rank loop. (Note that this definition 
is weaker than the one given in Subsection 14.21) . 

Theorem 15. A behavior P has a quantum representation of finite dimension d if and only 
if P admits, for some N , a certificate V with a rank loop, and rank(r 'n+ab) — d. 

Corollary 16. Let P be a behavior corresponding to a bipartite system where Alice's (Bob's) 
measurements have d& (ds) possible outcomes and such that each of the probabilities satisfies 
P(a,b) > 0, for all a £ A, b £ B. Let V be a certificate compatible with this behavior 
associated to the set of operators S, with S\ + ab Q S. Then, rank{T) = dAds implies that P 
has a quantum representation of dimension dAds- 

The proofs of Theorem [15] and Corollary [16] follow along the same lines as the proofs of 
Theorem [TO] and Corollary [TT1 

D Certificates and non-negativity of probabilities 

Let S\ + ab = {1, E a , E b , E a E b : a,b £ A, B} be the set of all sequences of length less than 1, 
together with all product operators consisting of one operator of Alice and one of Bob. The 
proposition here below states that the existence of a certificate T corresponding to a set of 
operators that contains Si + ab as a subset, thus in particular the existence of a certificate 
of order n with n > 2, implies that the elements P(a, b) of the behavior associated to T are 
proper probabilities, i.e., they are non- negative numbers. 

Before showing this, however, let us remind some notation introduced in Subsection 12.21 
We defined a behavior as a set of joint probabilities P = {P(a, b) : a £ A, b £ B} and 
implicitly assumed that they satisfy the no-signalling constraints P(a) = J2beY P{ a i b) an d 
P(b) = J2aex P( a ,b). To remove the redundancy associated with these constraints, we 
introduced the reduced outcome sets A and B so that P can be alternatively represented as 
P = {P(a), P(b), P{a, b) : a £ A, b £ B}. Having reminded this definition, it is now easy to 
see that the sets of operators S\ + ab = {1, E a , E b , E a E b : a, b £ A, B} and Sab = {E a E b : 
a £ A, b £ B} are linearly equivalent. 

Proposition 17. Consider a measurement scenario (A, B, X, y), and let P = {P(a,b) : 
a £ A,b £ B} be a set of real numbers. If there exists a certificate T for P corresponding to 
a set S such that Si+ab Q S, then the numbers P(a,b) represent proper probabilities, i.e., 
P(a, b) > 0, for all a and b. 

Proof. Let P admit a certificate as in Proposition [T71 Then, according to Lemma [7J P also 
admits a certificate associated to the set Si + ab, and thus also a certificate V associated to 
the set Sab = {E a E b : a £ A, b £ B}. Since T' y 0, its diagonal elements T' abab = P(a,b) 
must be non-negative. □ 



29 



E Proof of Lemma [131 



Proof. Suppose that there exists a pair of values (z,t), with \t\ < 1, such that 

M ^ =( nr o )h 0, (A-52) 



Q T R . 

with P = ^ ^ l^)' < ^ = (^x 1 x 2 )'^ = (t 1 ^ ecause 1*1 < ^ implies i? >- 0, 
Lemma [12] in section [5] states that the positivity of M 2jt is equivalent to the condition 
D = P — Q T Q~ l Q y 0. Now, D is a 2 x 2 matrix with non diagonal free entries and so its 
positive semidefmiteness is equivalent to demanding that D n ,D 2 2 > 0. Therefore, we can 
get rid of the variable z. Taking into account that t 2 — 1 < 0, we have that both conditions 
are equivalent to 

Pi<y< fa, (A-53) 



for a lj2 = X\%2 T a/ x i x 2 — x 2 — x l + 1j A, 2 — ^3^4 T a/ £3X4 — x\ — x\-\-\. 

It can be verified that |ai i2 |, l^i^l < 1. A solution max(«i,/5i) < t < min(a 2 ,/32) can be 
found if and only if 

«i < P2 

Pi < a 2 , (A-54) 

and the requirement that |t| < 1 translates into max(ai, min(o;2 ) P2) are n °t both equal 
to ±1. 

Now, it can be proven that condition (1A-54j) holds for any matrix M Z}t for which there 
exists a pair of values z, t that makes it positive semidefinite. To see this, notice that, for 
M Z)t to be positive semidefinite it is necessary that \t\ < 1. So, if such a couple of values 
exist, for any e > the matrix -^=(M 2ji + el) -7= is of the form (I3"gj) and there exists a 
pair of values {z! = z 1 j (1 + e), t' = t/(l + e)), with < 1, that make it positive semidefinite. 
Therefore, the vector (xi/(l + e) has to satisfy ( 1A-54I) . Because this holds for any e > 0, by 
continuity, also the vector (xj) will satisfy (1A-54H . 

Next we will prove that any vector satisfying (1A-54I) corresponds to a matrix of the type 
M for which there exists a couple of values (z, t) that make it positive semidefinite. Suppose, 
thus, that (1A-54j) holds. Two situations can arise: either max(«i,/?i) = min(a2,/52) = ±1 



or not. In the second case, we know that we can find a pair of values (z,t) such that 
M Zjt y 0, whereas in the first case it can be shown that x\ = x 2 = x; X3 = X4 = x' . 
But a positive semidefinite M matrix for this case is given by the formula M = D ■ 
(M* + diag(^ - 1, (^) 2 - 1, 0, 0)) ■ D, where D = diag(a;, x', 1, 1) and M* y is a 4 x 4 
matrix whose entries are all ones. Therefore, condition (IA-541) is necessary and sufficient to 
guarantee the existence of a pair of values (z,t) such that M z<t y 0. Making the change of 
variables x,- t — > sin(0j) in (IA-541) leaves us with fl3^|) . □ 
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